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INTRODUCTION 


The magnetic properties of Type II superconductors have been studied 
extensively by many researchers [i] 1 and, presently, a relatively good under- 
standing exists of the sources of hysteresis and magnetic instabilities in these 
materials. However, new effects continue to be observed 2 [2] and, in addition, 
neither theoretical nor empirical models have been found which can fully 
describe or explain more than limited amounts of published data [ 1] ., 3 

To contribute to a more complete understanding of the magnetic be- 
havior of Type n superconductors, an experiment was conducted to measure 
the way in which magnetic flux moves into and out of a sample under the 
influence of a changing applied field. Based on an extensive body of published 
information [ i] 4 and on prior experience [2] , we expected to obtain two types 
of flux motion information. First, as the external field H applied to such a 
sample is varied slowly, a relatively slow flow of flux through the sample 


1. See, for example, Reference 1 and the references therein. 

2. M. S. Lubell and E. W. Urban have independently observed large mag- 
netic flux flow oscillations apparently related to the onset of instabilities. 

3. See, for example. Reference 1 and the references therein. 

4. Ibid. 



wall occurs. This flux flow into or out of the superconductor is understood to 
produce spatial variations in the induction or flux density B within the mate- 
rial. Knowledge of flux flow characteristics would then permit inferences to 
be made of the distribution, motion, and causes of such flux density gradients. 

Second, as H changes over a wide range, the induced gradients 
spontaneously collapse at more or less regular field intervals, producing very 
rapid flux motion. During these flux jumps, as they are called, redistribution 
of flux occurs throughout much or all of the sample volume. Knowledge of 
the velocity of the jumps and conditions for jump initiation would give informa- 
tion concerning the thermal and magnetic diffusion properties of the super- 
conductor, and the occurrence of macroscopic flux instabilities. 

The work reported here concerns a series of flux flow measurements 
and the conclusions drawn from them. A schematic of the sample and sample 
holder arrangement used for the measurements is shown in Figure 1 and an 
overall schematic of the measurement instrumentation is given in Figure 2. A 
hollow cylindrical sample of Nb25%Zr with outer radius a and inner radius b 
was placed coaxially in the bore and at the midplane of a superconducting 
solenoid which provided the applied field H and permitted H to be swept at 
controlled rates or held constant. Hall effect probes near the sample mid- 
plane at the outer wall and in the sample bore directly measured H and the 
bore field H' respectively, following a technique used first by Kim and 
coworkers [3] . 


2 




I 



Fit 



2. Experimental apparatus. 








The rates of flux flow through the specimen walls, which we call 
4> = ^ , were measured by windings covering approximately three -fourths 
of the sample length. The single layer outer coil was wound directly on the 
specimen, while the multilayer inner coil was wound on a plastic insert shaft 
which also carried the inner Hall probe and three flux jump pickup coils. 
Descriptions of the high gain integrating digital voltmeters for measuring 
<p and 4> and of other measurements are given in Appendix A . Flux flow 
voltages typically ranged from a few to a few hundred microvolts, depending 
primarily on the field sweep rate H . 

Six flux jump pickup coils provided radial and axial flux jump motion 
information. Three identical coils were mounted near the outer wall, one at 
the midplane and one at each end; the other set of three coils referred to 
earlier were similarly placed in the tube bore. The six jump signals, typ- 
ically a few volts maximum amplitude and a few milliseconds duration, were 
displayed simultaneously on an oscilloscope and photographed. A represent- 
ative set of pulses for a flux jump is shown in Figure 3. It can be seen that 
the jump started near the bottom of the sample at the outside and proceeded 
upward and inward. 
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BACKGROUND DISCUSSION 


To provide a clearer understanding of the measurements to be 
described, a brief summary of the present understanding of the general 
magnetic properties of Type II materials is in order. Although all super- 
conductors are diamagnetic and strongly oppose the entry of magnetic flux 
into their interiors, Type II superconductors permit limited flux penetration 
in the form of flux quanta or fluxoids. The fluxoids can be visualized [4] as 
filamentary current vortices with small, normal phase cores which carry 

most of the field. The amount of flux in one fluxoid is 4> 0 = — — - ■ = 

-15 

2.07 x 10 weber, where h and e are Planck's constant and the elec- 
tronic charge, respectively. The field in a fluxoid decreases very rapidly 
with radial distance into the superconducting phase outside of the core, with 
the result that even at relatively high fields where fluxoids are densely packed, 
an appreciable portion of the volume of the material is superconducting and is 
available to carry nondissipative currents . 

If a sample of so-called ideal Type II material is placed in a slowly 
increasing field H, the following is understood to occur. At fields less than 
a lower critical field H , flux cannot enter the specimen because of the 
Meissner effect, a condition of perfect bulk diamagnetism. Above H , 

C 1 
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fluxoids nucleate at the sample boundary and move freely into the interior. 
Mutual fluxoid repulsion leads to a regular distribution of vortices throughout 
the ideal sample, their area density being the induction B. H is related to 
B by 

H = — B - M , ( i) 

Mo 

where M is the magnetic moment due to the spontaneous diamagnetic cur- 
rents which flow at the specimen surface. A plot of -M versus H for this 
situation is shown in Figure 4a. 

As H increases further, more of the specimen volume is occupied by 

normal fluxoid cores and -M decreases. Above the upper critical field 

H , the sample bulk is completely normal, and the magnetization is zero, 
cz 

Since flux flow is unrestricted, the curve of Figure 4a is reversible; an ideal 

Type II superconductor is not hysteretic. When H < H < H , the super- 

Cl cz 

conductor is said to be in the mixed state. 

In a real, disordered Type II material, such as the NbZr alloy used 
in this study, fluxoid motion is strongly inhibited by a process known as 
pinning [ 1] . Pinning centers have been identified as defects, dislocations, 
voids, precipitate particles, and many other types of micro structural features 
which are nearly unavoidable in technical materials. An important result of 
pinning is the existence of magnetic hysteresis, indicated in Figure 4b for three 
hypothetical materials. Pinning also leads to rapid spatial variations in flux 
density in such materials, as shown schematically on the right side of Figure 5 
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for the cylindrical specimen arrangement indicated in the inset. If, for exam- 
ple, the applied field has increased from zero to Hj, B varies in some man- 
ner from poH^ at r = a to zero at r = r*. Further increase of H can 
produce the succeeding induction distributions shown. 

Consider the Maxwell equation 


V x H - D = J 


( 2 ) 


For very slowly changing fields, D and D can be neglected. In our cylin- 
drical geometry, equation (2) then reduces to the scalar form 


cffi 

dr 


- Mo J 


(3) 


where J is a diamagnetic, azimuthal current density. Thus the local flux 
cffi 

density gradient is proportional to the local current density, and both 

vary with B. For example, when the applied field is Hj , J (B) , defined 

below, flows in the layer r t < r < a, as indicated in the left side of Figure 5. 

Two arguments demonstrate that as flux tries to overcome the p inn ing 

forces, it tends to move down the induction gradients. First, mutual fluxoid 

repulsion insures that motion will be toward regions of reduced flux density. 

1 dB 

Second, the Lorentz force F T = JB = - — - B — — existing between the local 

L p 0 dr 

current density and the individual fluxoids comprising B can be seen to be in 
the direction of decreasing B . If one calls the net local pinning force at a 
point in the superconductor Fp , then one may argue intuitively that flux will 
be almost completely pinned until F just exceeds F [5]. One can define 
a critical pinning condition as existing when 


il 



( 4 ) 



In this situation, maximum induction gradients exist and J(B) is called the 

critical current density J (B) . Then the magnitude of the Lorentz force is 

c 

J B and, if the functional behavior of F_ can be specified or inferred for a 
c P 

material, then the functional behavior of J (B) can be calculated, with the 

c 

aid of equation (4). 

On the basis of these last arguments, flux flow can be understood as 
the slow process of flux entry into or departure from a sample to establish 
the internal flux distributions, which are, in turn, governed by the critical 
current density function. Similarly, flux jumps result when the critical flux 
density gradients collapse. This leads to rapid redistribution of flux and to 
dissipation of the critical currents. For example, if a complete flux jump 
were to occur at H = H 2 of Figure 5, then at the conclusion of the jump, the 
field distribution would be as indicated by the dotted line, and would have 
vanished. Frequently in flux jump experiments it is found that the currents 
do not fully dissipate. Following such incomplete jumps, local regions of 
inhomogeneous flux distribution exist, together with the circulating currents 
which maintain them. 
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FLUX FLOW MEASUREMENTS 


Consider a flux flow measurement which starts from a condition of 

constant induction B. throughout the sample wall. The applied and bore 

fields are H. and H 1 . , respectively, and H. = — - — B. =H'. . Based on 

J J ] Mo J J 

the work of Coffey [6] and of Wipf and Lubell [7] , it was anticipated that 
this uniform field condition would be, in fact, well satisfied at the conclusion 

of most flux jumps because of the complete dissipation of the currents. As 

/ 

H is then increased or decreased by a small increment, a small amount of 
flux will enter or leave a thin annular region at the outer wall with some 
measured flow velocities ±<£j , respectively. If H continues to change 
at the same rate, then during the next instant flux must enter or leave a 
thicker annular layer and must fill or deplete the layer to a greater average 
induction increment than in the preceding moment. Therefore, the 
magnitude of the new measured flow rate ± will be greater than that of 
±$>! . Qualitatively, <fi will increase smoothly in magnitude with H, if 
the field sweep rate H is constant. The sign of 0 will be that of H , 
positive for increasing field and negative for decreasing field. The quanti- 
tative behavior of <p will depend on the exact nature of J(B) through 
equation (3) and its integrals. 
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A typical flux flow voltage variation measured for the NbZr sample 
with continuous, constant field sweep rate is shown in Figure 6. Flux flow 
voltage was recorded versus time, or equivalently versus H, since H was 
constant. Immediately after each flux jump, which appeared as a large 
voltage pulse, 5 the flux flow voltage dropped to a small value, increased 
rapidly, then rose linearly and much more slowly. There was no observable 
change in the form of the flux flow voltage, if the field sweep was stopped 
after each jump to insure that the sample was fully cooled to the helium bath 
temperature . 

The combination of the rapid flow voltage rise, followed by the slower 
increase, was not immediately understandable. A qualitative explanation of 
this behavior was proposed later and is described below. The interesting 
linear regions, each representing a constant flux flow rate increase , 
were found to vary systematically with field. In Figures 7 and 8 the slopes 
of these linear intervals are plotted as a function of the applied field at the 
midpoint of each interval for several values of H . Figure 7 shows <p for 
increasing H . The flux flow slopes have a maximum at about 6 kilogauss, 
dip through a flat minimum, then rise smoothly. For decreasing H (Fig. 8) , 
no peak appears. 


5. Note that this integrated pulse is simultaneous with, but obscures the 
details of, the type of pulse data shown in Figure 3. 
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Figure 6. Typical flux flow 


H = 9.8 kG 



tage data; continuous field sweep mode. 









FLUX FLOW CALCULATIONS 


Following the observation of these regular flux flow rate features, 
we sought to verify them analytically. A calculation was made of the flux 
distributions and flux flows to be expected in a cylindrical specimen. It was 
assumed that each smooth flux flow sequence would start from a uniform field 
condition = p 0 H. , as in Figure 9. It was further assumed that when H had 
changed by some amount H (less than that required to cause the next jump) , 

3 . 

one of the two quasistatic situations shown schematically in Figure 9 would 
apply with 

H = H. ± H , (5) 

J a 

the upper and lower signs referring to H > 0 and H < 0, respectively. 

r . is the radial distance to which the change in B(r) has penetrated; it 
mm 

depends on H. and H , and on the particular J (B) . 

j a c 

In order to guarantee that the induction distribution B(r) was indeed 
quasistatic, one had to assume that the processes of thermally activated flux 
creep, described originally by Anderson [8] and Kim et al. [3, 9] , could be 
neglected. Thus if the field sweep were stopped at any level, the B(r) at 
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Figure 9. Induction distribution and calculation notation. 









that instant would cease changing, and $ should very quickly become zero. 
Restart of the sweep at the same rate should then cause <£ to return imme- 
diately to the same amplitude and slope which existed when the sweep was 
stopped. On the other hand, the occurrence of appreciable flux creep after 
the sweep was stopped would cause a gradual (j> decrease to zero and would 
allow B(r) to change somewhat, giving rise to a different amplitude and 
slope upon sweep restart. During the series of experiments, it was verified 
several times that a sweep stop, followed by several seconds wait, then 
sweep restart did, in fact, give the behavior expected for a quasistatic situa- 
tion in which flux creep was negligible. 

To indicate the complete functional dependencies of our problem, we 
rewrote equation (3) in the form 


dB(r, t) 
dr 


- Mo J c [b ( r , t)] 


( 6 ) 


It remained only to specify a functional behavior of and to integrate 

equation (6) to obtain the variation of induction with radial position in the 

specimen. The boundary conditions are B(a) = (H. ±H ) and 

J a 

B(r . ) = (UqH.. 
min u j 

Currently eight essentially different models for the field dependence 
of critical current density of Type II superconductors have been proposed in 
the literature. They are summarized in Table i, with references to the 
works in which they appear. In each model, ol is assumed to be a 
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TABLE i. CRITICAL CURRENT DENSITY MODELS 


Model 

J 

c 

References 

Bean- London (BL) 

a 

10 , 11 

Kim (K) 

a 

B 0 + B 

12 

Anderson-Friedel-Silcox (AFS) 

a 

b" 

8, 13, 14 

Yasukochi (Y) 

a 

B 1 / 2 

15 

Irie (I) 

a 

B 7 - 1 

16 

Fietz ( F) 

“ exp (‘t) +Y 

17 

Goedemoed (G) 

“* M 0 H c2 ' B ) 

18 

Alden-Campbell-Coffey (ACC) 

"o H C2 : B 
a 

B 1 / 2 

19,20, 21 


temperature dependent constant. 6 For the Kim model, B 0 is a constant 
whose physical meaning is not clear. Empirically it has been found that B 0 


6. For the following reason, temperature effects were not considered in the 

calculations. Experimentally, the sample was fairly completely exposed 
to the liquid helium bath, so that if the thermal conductivity of the super- 
conductor were large enough and the rate of heat generation due to dis- 
sipative flux motion were small enough, a constant sample temperature 
equal to the bath temperature 4. 2®K could be expected. Wipf and Lubell 
[7] state that for a Nb25%Zr sample approximately of the present 
specimen size, a uniform, constant temperature condition will exist 
throughout, so long as H is less than about 10 4 G/s . Since maximum 
H was 200 G/s , the isothermal assumption was justified. 
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is of the order of a few kilogauss for NbZr and NbTi alloys and Nb 3 Sn. 

For the Irie model, y is a number which should lie in the range i^y^2. 

Having obtained from the integration of equation (6) an expression for 
B(r,t) or, as is the case for some of the models, an expression for r(B,t) , 
the instantaneous total flux (p (t) in the sample can be calculated. Since 
induction B is flux density or flux per unit area, the total flux in the speci- 
men wall is given by an integral of B over the area of the cylinder cross 
section, 


a 

0(t) = 27 r f B(r,t) rdr 
0 

r . 

min a 

= 2 7t f n o H. rdr +27 r f B(r, t) rdr 

0 J r . 

mm 

a 

0(t) = 7T JU 0 H. r 2 . + 27T f B(r,t) rdr 

™ l mm J 

J r . 

mm 


(7) 


(8) 


The first term in equation (8) represents the total flux within the cylindrical 
region 0 < r < r . due to the field H. of the previous jump. Since r 


mm 


J 


mm 


decreases with increasing H , this term is a function of time. The second 

term is the total flux in the annular region r . <r^a in which the 

mm 

induction has changed from H . 

If the specimen is a solid cylinder and b = 0 , equation (8) is valid 

only if r (H ) > 0 . For greater H , flux reaches the axis from all 
mm a a 

sides and a reduced flux flow rate is required to establish a given induction 
on the axis. 
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If the specimen is hollow, equation (8) is valid only if r > b . For 

mm 

greater H & , flux must flow into and across the sample wall at an increased 
rate to permit the bore field to increase smoothly. 

Our primary concern in these calculations was the determination of 
0( t) and 0(t) . It would, therefore, have been possible to avoid a major 
calculation step by differentiating equation (8) directly. However, since the total 
flux in the sample as a function of field could be found from the experiments 
simply by measuring the area under the flux flow rate curves, an explicit 
expression for 0(t) was considered valuable. Consequently, except for the 
Fietz (F) model, 0( t) was determined. It is clear that 


0(t) = 


90 dir 
9H dt 


90 


9H 


H 


(9) 


and that 


•• 90 dH 90 

* (t> = m dT + IT 


( 10 ) 


• 0 (b 0 (p • • 

However it is assumed that H is constant so that — = — rr H vanishes. 

dt dH 


By means of equation (9) , it is seen that equation (10) reduces to 


0(t) = 


9T? 


H^ 


( 11 ) 


In our calculations we defined the following parameters: 

H. uniform field in the sample after each jump, 


H (t) 
a 


units are ( — — kG 
\ Mo 

applied field change from H. 
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H(t) = H.± H (t) 
J a 

H> ± H 

a 


B(r,t) 
a = + a. 

c 

<p(t) 

<Pi t) 

0( t) 


x = 



H. ± H 



x. = 
3 





B 


Mo 


H 


c2 


total applied field kG^ 

rate (and direction) of applied field sweep 



flux density or induction in specimen (kG) 

proportionality constant (dimensions given 
for each model) 

total flux in specimen (weber, Wb) 

rate of change of flux or flux flow rate 
(Wb/s) 

flux flow acceleration or flux flow rate slope 
(Wb/s 2 ) 

dimensionless applied field 
dimensionless jump field 
dimensionless induction 


mm 


sample outer radius (meter) 

minimum radius at which B has changed 
(meter) . 


The upper and lower signs in H and H account for the increase or 
decrease in field for positive and negative field sweeps, while the signs in 
a similarly account for the diamagnetic nature of the superconducting cur- 
rents. In the tables, J £ , r^ , B(r,t) or r(B,t), (pit), < pit ), and 
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(S(t) are shown. It should be understood that r . =r . (t) through the time 
^ min min 

dependence of H , and that the three flux functions depend on H. and H , as 
a j a 

well as time. 

In Appendix B are listed the results of the flux flow calculations based 
on the eight critical current density models of Table 1. In each expression 
for cp(t) the first term is seen to be the total flux contained within the outer 
wall due to H.; for each inter jump interval this term is independent of 
time. 

Although calculated flux flow rates will not be shown explicitly here, 
it may be stated that the <p results of all of the eight models predict a smooth 
rise of the flux flow rate from zero following each flux jump. There is no 
suggestion that the sharp breaks in the curves of Figure 5 should have been 
observed. This clear contradiction between the calculations and the measured 
flow rates led one to question the assumptions made in the calculation. Wipf 
suggested that the problem might stem from the acceptance of a condition of 
field uniformity in the sample following each jump. If, in fact, large 
inhomogeneities remained, then the flux would not necessarily enter smoothly 
nor assume a well-behaved distribution. 
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STOP-HEAT FLUX FLOW MEASUREMENTS 


Since it was not realistic to try to modify the calculations to account 
for some unknown flux distributions, the experiment was modified to insure 
that the constant induction assumed to exist at the start of each inter jump 
region would indeed exist. A single layer, nichrome wire heater was wound 
noninductively over the outer flux pickup coil and on the central pickup coil 
shaft. It was then possible to heat the entire specimen to above its transition 
temperature at each jump field to permit the field throughout the wall to 
become uniform and equal to the applied field EL . 

The measurements were then made in the following manner. Immedi- 
ately after each jump the field sweep was stopped. The sample was then 
heated for approximately 5 seconds at a heater power of approximately 10 
watts, which was sufficient to allow the flux to reach a quiescent, redistrib- 
uted condition. Finally, the specimen was allowed to cool for 10 seconds to 
insure that its temperature had returned to 4.2°K, then the field sweep was 
restarted. Typical data for this stop-heat mode of operation over one inter- 
jump interval are shown in Figure 10. There is a smooth rise of the flow rate 
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H = 186 


sec 




>; 1 division = .80 sec = 149 G 

;e data; stop-heat mode. 


from zero as predicted by the calculations . The flow rate slope 0 is 
approximately constant. 7 These results confirm the supposition made in the 
preceding section that, following each jump, the unheated specimen contains 
induction inhomogeneities. 

Figures 11, 12, and 13 show flux flow rate slope data for the stop- 
heat measurements. The average slope in weber/ sec 2 of each approximately 
linear-jump flow voltage curve is plotted versus H , the applied field at 

3 , 

the midpoint of the interval. Figure 11 shows the discrete data for three 
sweep rates and both sweep directions. Figure 12 gives the H =49 G/s 
points of Figure 11 on an expanded scale. Figure 13 is a semilogarithmic 
plot in which curves summarize the data of Figure 11. The slopes are 
represented by a smoothly increasing function of average applied field (the 
low field peaks of Figure 7 do not occur) . 


7. At fields above l5kG (approximately) , the data have a moderate downward 
curvature; in all cases the flow rate slopes are taken from linear inter- 
polations between the endpoints of the smooth regions. 
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Figure 11. Flux flow rate slopes; stop-heat mode; 
all sweep rates. 







COMPARISON OF CALCULATIONS AND DATA 


In order to compare the 0 formulas of Appendix B for the several 
critical current density models with the experimental results, the following 
procedure was used. Actual experimentally measured fields were used in 
solving the 0 equations: H was taken to be the applied field observed 
following each heating and recooling; H was one-half of the subsequent 

3 , 

constant slope interjump interval. Thus H = H. ± H was the midpoint 

J a 

field of each interjump region. The constant a was chosen for each model 

so that the calculation fit the H = 118. 5 G/s curve at point A of Figure 12. 

This fit point was chosen, since it applies to an intermediate field, where 

presumably all of the models should be reasonably applicable, and it applies 

to an intermediate value of H . For point A the parameters have the 

following experimental values: 0=0. 0971 Wb/s 2 ; H = 118. 5 G/s; 

H. =20. lOkG; H =1.53kG; H =21.63 kG. 

J a 

The results of the 0 calculations are shown in Figures 14 and 15. 
The same 118. 5 G/s data have been replotted three times in each figure to 
separate the various curves for clarity. Curves for three values of B 0 for 
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Figure 14. Flux flow rate slope calculations. 





Figure 15. Flux flow rate slope calculations. 



the Kim model are shown (B 0 = 0 gives the Anderson -Friedel-Silcox (AFS) 

result) . For the Irie calculation, which reduces to the Bean-London, 

Yasukochi, and AFS curves for y = 1, 3/2, and 2 respectively, only the curves 

for y = 2. 5 and 3 are shown. Note, however, that from Irie's original 

work there is no immediately apparent physical explanation for y > 2. For 

the Goedemoed and the Alden-Campbell-Coffey (ACC) calculations, H „ was 

c2 

taken as 63.6 kG, a measured value discussed in Appendix A. 

The <p computation for the Fietz model with its three adjustable 

parameters was made as follows. It was required that J c go to zero at H ^ , 

a necessary condition which Fietz et al. did not use in their work. This 

/ ^ H c 2 \ 

gave the result that y = - a exp I - — ) , leaving two parameters oi 

and /3. In their paper Fietz et al. empirically found (3 to be of the order of 
5 kG for two coils wound of Nb25%Zr wire. Therefore, (3 = 5 kG was chosen 
and the computation was fit at point A as usual, giving the strange behavior of 
Figure 15. Larger values of )3 were then chosen, in order to obtain a fit to 
the data at point B of Figure 12, as well as at point A. The very good fit 
in Figure 15 resulted when (3 = 25 kG . This value is so far from those 
measured by Fietz that the agreement must be considered as physically 
inexplicable as is the agreement of the Irie computation for y > 2 . 

The apparent variation of the Bean-London (BL) curve with H is 
spurious. The BL expression for depends on H , but not on Eh, so 
the experimental variation of the size of the jump intervals leads to an 
apparent variation of <£> . This effect undoubtedly affects the curves of the 
other models but to a much smaller degree. 
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If the Fietz model is disregarded because of its unrealistic value of 
(3 , ACC results in the best empirical fit over the widest field range. Never- 
theless, the calculated <p is 30 percent high at the lowest field jump and 
about 8 percent low at the highest field jump. 
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NEW CRITICAL CURRENT DENSITY MODEL 


The unsatisfactory agreement with the present experiment of the 

flux flow calculations for all of the eight published critical current density 

functions (Table 1) led us to seek a better functional form for J . Since 

c 

it appeared that a probable source of the failure of the models might be the 
assumption that a is independent of field, the variation of a with B was 
determined explicitly. Choosing the Kim model as a starting point, the 
successive values of a required to match the <j> formula to the series of 
measured interjump flow voltage variations were calculated. These ads 
are plotted against average applied field in Figure 16. By means of a least 
squares curve fit they are shown to lie approximately on a straight line which 
approaches zero at about 61 kilogauss. Because this maximum field is of 
the order of magnitude of H , it appeared reasonable to assume that 
o'(B) for the Kim model varies as 



where a is a field independent, temperature dependent constant. Inserting 
c 

equation (12) into the Kim function of Table 1 gives 


J = a 
c ' 


H - B 
c2 

c B 0 + B 


a functional form which is seen to vanish at the upper critical field H and 


to approach a large, finite value a 


^0 H c2 


e B n 


at B = 0. 


The proof of the validity of this model lay, of course, in its ability to 
explain the present flux flow measurements and to agree at least reasonably 
well with other published results. Therefore, the total flux, flux flow, and 
flux flow slope calculations were made from equations (6), (8), (9), and 
(13). The results are summarized in Appendix C. H was given its 

C £i 

measured value of 63. 6 kG (Appendix A) . Here, we have two adjustable 

parameters, a and B 0 , which were chosen to give exact empirical agree- 
c 

ment of the </> calculation with points A and B of Figure 11; for this fit, 
o? c = 3. 3 x 10 3 A/cm 2 and B 0 = 6. 5 kG. (p curves were then calculated 
for all of the flux flow slope data of Figures 11, 12, and 13. The results 
shown in Figures 17 and 18 agree remarkably well with the measurements. 
We may, therefore, conclude that equation (13) is a very good representa- 
tion of the field dependence of the critical current density, at least for 
Nb25%Zr. It remains, of course, to determine the validity of equation (13) 
for other Type II materials. 

It is possible to verify the choices of <x and B 0 made from the flow 
slope considerations by means of an independent experimental measurement. 


In Figure 5 it is observed that when H = H , the flux distribution B (r) just 

s S 
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reaches the tube bore . It follows that r . =b at that field. Subsequent 

mm 

increase of H causes flux to enter the bore and H' begins to increase. 

During the flux flow experiments, a similar sequence occurs each time the 

flux front starts from Eh and reaches the tube bore without a flux jump 

having previously disturbed the B distribution. Figure 19 shows an 

experimental H' versus H curve taken for H = 20 G/s. At points E, G, 

and I, the flux front clearly reached the bore and H' increased, if only 

briefly. From the expression for r(B) in Appendix C and using the original 

values of H , a , and B 0 , we calculated the radial variation of B for 
cz c 

selected values of H and H to match the fields on the increasing field 

j a 

portion of Figure 19. The B results are shown in Figure 20. Note that the 

calculated B curves reach the bore at the exact fields observed in the 

experiment. This confirms the empirical choice of a and B 0 from the 

c 

fit to the flux flow slope results. Note also that flux jump number 2 initiates 
just as the flux front reaches the bore wall. This last observation directly 
confirms earlier suggestions [5] that a hollow superconductor is especially 
susceptible to flux jumps just when a sudden increase of flux flow must occur 
to permit the entire bore to begin to fill with additional field. Heating in the 
specimen, resulting from the more rapid flow, reduces the pinning forces 
and a jump may ensue. 
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Figure 19. X-Y record of bore and applied fields; 
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Figure 20. Induction in tube wall for experiment of Figure 19; calculated from new model. 



Thus far, solutions to the very complicated formulas for cp 

given in Appendices B and C have been obtained. It is instructive to consider 

the behavior of the critical current functions which served as the starting 

point of each flux flow calculation. Using the values of the various parameters 

which gave the curves of Figures 14 and 15, the J curves for all of the 

c 

models except that of Irie are plotted in Figures 21 and 23. If J for the 

c 

new model is taken as the correct value for the present sample, the deviations 

of the different J models can be seen to reflect the deviations from the 
c 

experimental results in Figures 14 and 15. 
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Figure 22. Critical current density functions with proper 
high field behavior. 


4 ' 


DISCUSSION OF THE NEW MODEL 

Equation (13) was found empirically and a physical basis for it remains 
to be established. In particular the significance of the constant B 0 is as 
obscure for the new model as it has been since it was originally introduced 
by Kim [9] (Table 1) . It is striking, though perhaps only fortuitous, that 
the value of B 0 is just that field at which the flux flow rates for the con- 
tinuous field sweep have a relative maximum (Fig. 7) . 

It is interesting to compare the Goedemoed and Alden -Campbell -Coffey 
critical current formulas with the new one (N) . At low fields where 

<B«jn 0 H the three current density models vary as jUoH , 

C 1 CZ cz 

-i/2 -i 

B , and /jqH (B 0 + B) for G, ACC, and N respectively. Above 

CZ cz 

H . , which is of the order of a few hundred gauss for Nb25%Zr, the 
cl 

Lorentz force or pinning force functions then vary as /u 0 H B, £< 0 H B 1 ^ 2 , 

c z cz 

and u,H ( ) respectively. The low field importance of B 0 in the 

c2 \ Bq + B / 

new model is apparent. 

At high fields where B 0 may be neglected, the three forms may be 

summarized as J ~ (/i 0 H -B)B n , where n=0, l/2, 1 for G, ACC, 
c cz 

and N respectively. 
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The n dependence can be seen in Figure 22 near H . The Lorentz force and 
the net pinning force at high fields are then 

F L =F p~ ( ^° H c2- B)Bl " n ; B>>B o > U 4 ) 

and for N they vary as juqH - B. Coffey [21] has calculated that the 
fraction of the specimen cross section which is outside the fluxoid cores and 
remains superconducting is also proportional to ju 0 H - B. Consequently 

C M 

only the new model is seen to be consistent with the quite reasonable argument 
that the net pinning force on the fluxoids varies directly with the local area 
density of superconducting phase in the material. 
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SUMMARY 


In this report we have discussed new observations of magnetic 

flux flow behavior in the Type II superconductor Nb25%Zr. The rate of 

flux flow <^> into or out of a cylindrical specimen increases linearly with 

field in the interval between flux jumps, provided the process in each interval 

starts from a condition of constant induction throughout the specimen. The 

slopes of the linearly increasing flow rates, which represent $, are 

themselves increasing functions of average applied field. Since the flux in 

the sample and its derivatives with respect to time can be calculated from 

an assumed variation of critical current density J (B) , it is possible to study 

c 

the applicability of various functional forms of which have been proposed 

in the literature. Although some of the eight previously advanced models 

were able to explain the present data only over limited ranges of applied 

field, a new empirical model proposed and discussed in this report is 

shown to agree with the experimental results at all fields, except near H 

Vo H c2 _B C 

where anomalous behavior occurs. The new model, J = a — — — — — — , 

c c B 0 + B 

is well behaved at very low fields and vanishes properly at H At 

C Zj 

intermediate fields, its variation closely resembles short sample critical 


50 



current densities reported in the literature. Although the physical signifi- 
cance of the constant B 0 remains obscure, its empirical value 6.5 kG is 
close to the values found for this parameter in the simpler, related Kim 
model. It remains to provide an explanation of the new model from considera- 
tions of microscopic flux pinning interactions in the superconductor. 8 


8. Portions of the work discussed here are also reported in Reference 22. 


M438 
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APPENDIX A 

EXPERIMENTAL DETAILS AND ADDITIONAL MEASUREMENTS 


Because some relatively novel experimental problems were encoun- 
tered in this study and because the special dc amplification system used in 
the flow measurements has very good accuracy, sensitivity, and general 
utility, a more detailed discussion of the apparatus will be given. 

The flux flow measurement system employed for this study was in 
effect a dual, rapid, high gain integrating digital voltmeter with simultaneous 
digital and analog output records . One channel was used for the outer flux 
flow pickup coil, the other for the inner coil; one channel is described. All 
of the apparatus in this dc amplifier was standard Hewlett-Packard instru- 
mentation, but its use in the present system combination is noteworthy. A 
2212A Voltage-to-Frequency Converter (VFC) converted the low level dc 
signal into a pulse train whose repetition frequency was directly proportional 
to the input voltage. The VFC sensitivity was 10 Hz per pV. The VFC 
pulse frequency output was counted by a 532 IB Frequency Counter at a rate 
of 10 frequency readings per second. These frequencies corresponded to the 
average input voltage in microvolts during the counting time. They were 
digitized and recorded on one of two channels of a 5050B Digital Printer. A 
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temporary data storage feature in the printer reduced the time during which 
the printer inhibited the counter to 100 ps per print cycle. Consequently 
the system monitored the input signal 99. 90 percent of the time and the 
probability was very small that any significant data fluctuations, including 
flux jumps, would be missed. 

As they were recorded digitally, the frequency data were also fed to 
a 581A Digital-to-Analog Converter (DAC) which generated a dc voltage pro- 
portional to the number represented by the three least significant digits of 
each counter reading. This voltage was recorded on a 7100A Strip Recorder 
whose full scale deflection was adjusted to match the DAC output voltage 
range for a 0 to 999 reading. Thus, although the counter readings were often 
greater than 999, the recorder pen always remained on scale by recycling 
as often as necessary. Note, for example, the 636 pV and 1778 pV deflec- 
tions in Figure 10. It can be seen that this method effectively expanded the 
width of the final analog recording, while maintaining the same absolute 
sensitivity. Ambiguities introduced by pen recycling were removed later by 
reference to the digital record. Introduction of the 0. 1 second integration 
time into the data produced small, unobtrusive steps in the final analog re- 
cord, as was seen in Figures 6 and 10. 

To eliminate thermoelectric voltage problems in the flux flow 
measurements, the pickup coils were formed of a continuous conductor 
from one room temperature VFC input through the coil winding in the liquid 
helium and back to the other VFC input. The VFC-to-coil connections 
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themselves were enclosed in a cardboard box to protect them from rapid 
temperature changes . Drifts in the entire system were very slow and could 
easily be corrected for by occasionally rezeroing the VFC's. 

Although the field sweep rates were automatically and very accurately 
controlled at intermediate magnetic fields by the magnet power supply, 
problems were encountered at low (<8 kG) and very high (>36 kG) fields. 
In these regions for H > 0, the magnet sweep rate changed in a slow, 
reproducible manner because of changing magnetic inductance. These 
anomalous sweep rates were independently measured by means of a bare 
copper coil in the solenoid bore and the data plotted in Figures 11 through 14 
and 16 through 17 were corrected accordingly. The data in these figures 
accurately reflect the flow slope behavior which would be observed for the 
stated constant H values. 

Since one parameter of three of the important critical current models 

considered in this report is the upper bulk critical field H , its direct 

c2 

measurement was necessary. Published values [23] of H for Nb25%Zr 

C 

are spread around 70 kG and have large uncertainties. Since H is a 
function of composition, and the exact composition for this specimen was not 
known, an independent critical field determination was made. Under normal 
circumstances H is the field at which the sample becomes normal and the 

C zj 

difference H' - H between the bore and applied fields becomes zero. Since 
both branches of the critical state curve converge to H' = H at H , the 
upper critical field is observed as the field at which the two branches merge 
with the H 1 =H Hall probe characteristic. 
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The NbZr sample studied exhibits an interesting property known as 


the peak effect. It is common to many Type n materials and has been attrib- 
uted [24] to increasing resistance to mutual interpenetration by the flux lines 

in the sample near H „ . Its existence necessitates the definition of two 

c2 

upper critical fields. High field H 1 versus H data are shown in Figure A-l 

in which the peak effect is quite evident; at several fields on both critical 

state curves the slow sweep was stopped and the sample heated to above T to 

c 

define the Hall probe characteristic. 

Above approximately 55 kG for H > 0 the field difference H' - H, 

instead of continuing to converge monotonically toward zero, begins to 

increase with field. It passes through a maximum and then falls rather 

rapidly to zero. For H < 0, a similar maximum in H’ - H is seen. This 

rapid high field convergence of the critical state envelopes permits relatively 

precise determination of the true bulk upper critical field, which will be 

simply referred to as H One would expect that any reasonably simple 

cz 

physical model which explains the magnetization and flux flow data at fields 
below 55 kG should extrapolate smoothly into the region in which the peak 
effect is observed. It follows that the value of upper critical field which must 
be used in such a model is that field toward which the magnetization is 
converging when the peak effect commences. This effective bulk upper 
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as seen in Figure A-2. The peak effect is very pronounced and at its high 
field side drops steeply and smoothly to an intercept at 69.4 ± 0. i kG. A 
very linear decrease occurs prior to the onset of the peak effect rise. A 


linear least squares curve fit leads to an extrapolated effective value 


H c2(eff) 


of 63. 6 ± 0. 5 kG, the value which was used above in the fit of the critical 


current density models to the flux flow data. 
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APPENDIX B 

CALCULATION RESULTS - PUBLISHED MODELS 
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Goedemoed Model 
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APPENDIX C 
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